Microscopic Theory of Damon-Eshbach Modes in Ferromagnetic 

Films 



R. V. Leite and R. N. Costa Filho 
Departamento de Fisica, Universidade Federal do Ceard, Caixa Postal 6030 
Campus do Pici, 60451-970 Fortaleza, Ceard, Brazil 
E-mail: rai@fisica.ufc.br, FAX:-l-55-85-288 9903 

Abstract 

The surface spin wave branches in ferromagnetic films are studied using a 
microscopic theory which considers both magnetic dipole-dipole and Heisen- 
berg exchange interactions. The dipole terms are expressed in a Hamiltonian 
formahsm, and the dipole sums are calculated in a rapidly convergent form. 
The Damon-Eshbach surface modes are analyzed for different directions of 
the spin-wave propagation and also for different ratios of the strength of the 
dipole interactions relative to the exchange interactions. Numerical results 
are presented using parameters for Fe and GdCls. 
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It is well known theoretically that, in quasi-two-dimensional systems, such as ultrathin 
ferromagnetic films, the short-range exchange interactions alone are not necessarily sufficient 
to establish a Ferromagnetic ordering above the ground state, and that it is necessary to 
take into account the anisotropy and long-range character of the dipolar interactions [|I|. 
Regarding the spin dynamics in ferromagnetic films, a macroscopic (or continuous medium) 
theory for the dipole-dominated regime was given by Damon and Eshbach in terms of 
magnetostatic modes. They identified a surface branch of the spectrum, now known as the 
Damon-Eshbach (DE) mode. Recently, microscopic theories have been used to study the 
dipole-exchange spin waves (SW) in ferromagnetic and antiferromagnetic films 

The DE mode has a number of striking and unusual properties. For example, the exis- 
tence condition of this mode depends on the propagation angle relative to the principal axes 
of the crystal. The x-direction in our case, being within the range 6 < 6^, where 6(. 
depends on the characteristics of the system, such as the applied field and dipolar strength 
(see Ref. 6 for a review). The aim of this paper is to study the DE modes in ferromag- 
netic films at low temperatures (T <^ T^, where T^. is the Curie temperature) by using a 
microscopic theory, and including both magnetic dipole-dipole and Heisenberg exchange in- 
teractions. The dipole terms are expressed in a Hamiltonian formalism, and the dipole sums 
are calculated in a rapidly convergent form. In order to find the critical 6c where the DE 
modes vanish, we have to study the dependence of the DE modes and the SW energies (for 
each of the discrete branches) on the variation of the exchange parameter, on the strength 
of the dipolar interaction, and on the wave vector direction of propagation. Another impor- 
tant factor in the analisys is the effect of different dipole interaction strengths (relative to 
the exchange)on the dispersion relation of the SW. Numerical results are presented using 
parameters for Fe and GdCla. 

Consider a ferromagnetic film with atomic layers arranged on a simple cubic lattice, 
with lattice constant d. The wave vector q makes an angle 6 with respect to the x-axis, and 
the Zeeman field Hq is assumed to be parallel to z direction. The surface of the film is in the 
X — z plane, while the y-axis is perpendicular to the film. The Hamiltonian can be written 
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as 

^ = E J^A ■ - gf^BHo E 5f + ligf^B)' E D^fs^s^, (i) 

ij i ij 

where Sj is the spin operator at site i and Jij is the exchange couphng between sites labeled 
i and j. We assume that the exchange coupling is Ji and J2 for nearest and next-nearest 
neighbors respectively. The second term of the above equation is the Zeeman term, and 
the dipole-dipole interaction between spins is represented by the last term, where a and f3 
denote the components x, y, or z, and 

< = {|r.,f'5<./3-3r^rg}/|r.,f (2) 

whit Tij = Tj — Fj, and the case i = j excluded from the sums in Eq. (1). 

Using the Holstein-PrimakofF transformations, the Hamiltonian can be expanded as if = 
ff(2) _|_ ^(3) _|_ ^(4) _|_ _ (apart from a constant), where if^™) denotes the term involving 
a product of m boson operators. Here we are considering the non-interacting (linear) SW 
modes that are described by the quadratic term H^'^\ Therefore, the Hamiltonian takes the 
form: 

H = ^{Ann'{q)al^agn' + (q) [ag„a_5„' + aJ„aLg„/]}. (3) 

qnn' 

Here we are using a representation of the boson operators and a in terms of a 2D wave 
vector q = (g cos 9, q sin 9), q being the modulus of q, parallel to the film surface and indices 
n, n' (= 1,2, ...,A^) that label the atomic layers parallel to the surface. The coefficients 
^nn'i^l) and i?„„/(q) are described in Ref.5. A linear transformation to new boson operators 
can then be found such that H takes a diagonalized form 0. Its eingenvalues yield the 
energies -Ei/(q) of the discrete SW branches, where i/(=l,2,...,N) labels the branches. 

In Fig. 1(a), we show the behavior of the lowest SW branches for very small wave 
vectors in a GdCls film with = 20. For this material, we use the parameters; AnM = 
{gfiBS)/d^ = 0.82 T to characterize the dipolar strength, the bulk exchange field Hex = 
6S{Ji + 2J2) = 0.54 T, and Hq = 0.36 T. We have taken the case of ^ = (the Voigt 
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geometry), for which it is known that there is a Damon-Eshbach (DE) surface mode 
in the magnetostatic continuum hmit where exchange is neglected. It can be seen that a 
purely magnetostatic DE mode starts at = g = at the frequency 18.25 GHz. There are 
hybridizations until the eighth mode, and the DE mode tends to be flat around the value 
21.5 GHz with increasing q^. On the other hand, in a material such as Fe that has a much 
smaller value for the ratio AttM/Hex, it is necessary to consider a much larger film thickness 
in order for the DE mode to hybridize with the lowest bulk branch. This is because in a 
thicker film, the density of modes increases and create more modes in the region where the 
DE mode propagates, making hybridizations possible. An example for Fe is given in Fig. 
1(b) considering = 120, AtiM = 2.14 T, Hex = 2140 T, and Hq = 0.054 T. Only the 
lowest branches of the spectrum are shown in Fig. 1(b), where it can be seen that much 
smaller values of Qx are necessary (compared with the previous example) to see the analogue 
of the DE mode in the spectrum. 

In Fig. 2, we show results for the SW frequencies of the above-mentioned materials 
considering the SW propagation direction 9. In Fig. 2(a), the dependence of the lowest 
GdCls SW branches on 6 is shown for a fixed value qa/n. For small wavevectors, the 
branches decrease their frequency values as ^ = goes from to ^ = 7r/2. As in Fig. 1, 
the hybridization of the modes is seen to decrease until the second mode. The two lowest 
branches are the most affected by the SW propagation direction. The others tend to become 
fiat after certain angle values. The figure also shows the reciprocity of the dispersion relation 
as the angle goes to vr. The same effect can be observed in Fig. 2(b) considering Fe. In this 
case, the dependence of the modes on the angle is much weaker than in the case of GdCla, 
and only the two lowest modes are affected by the propagation angle. 

In conclusion, we have used a microscopic theory to study Damon-Eshbach SW modes 
in the dipole-exchange regime in thin ferromagnetic films. We analyzed the dependence of 
the DE modes on the SW propagation direction and the ratio between dipolar and exchange 
strengths. The results show that the DE modes tends to vanish as the system change from 
the Voigt geometry to a general propagation direction, and that this change depends on the 
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strength of the ratio between the dipole and exchange interactions. 
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FIGURE CAPTIONS 



Fig.l The linear SW dispersion relation for the lowest branches of (a) a 20- layer GdCla 
film and (b) a 120-layer Fe film. In both cases ^ = 0. 

Fig. 2 The SW frequencies as a function of the propagation direction 9 for (a) a 20-layer 
GdCla film taking qa/n = 0.02 and (b) a 120-layer Fe film for qa/n = 0.002. 
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FIGURES 




FIG. 1. R.V. Leite at al. 
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